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MODEL  STUDIES  OP  SEISMIC  ENERGY  DISTRIBUTION 
AROUND  DIFFERENT  TYPES  OF  SOURCE 

Progress  Report  No.  4 
for  the  quarter  ending  November  30,  1961 

I.  SIMULATED  MODEL 

One  of  the  objectives  of  this  research  has  been  to 
accumulate  a  "library"  of  pictures  showing  the  patterns 
of  first  motion  around  different  types  of  source  as  a 
function  of  azimuth  for  a  uniform  medium.  The  gathering 
of  data  for  this  purpose  has  now  been  completed.  Data 
on  first  shear  (S)  waves  under  the  same  conditions  has 
also  been  recorded.  Four  source  systems  simulating  an 
explosion,  a  simple  unidirectional  force,  a  single 
couple  and  a  double  couple  have  been  modeled. 

The  patterns  are  modified  in  the  presence  of 
velocity  discontinuities  near  the  source.  Five  typical 
conditions  have  been  studied:  a  boundary  passing  thru 
the  focus;  one  to  one  side  of  the  focus;  two  boundaries 
symmetrically  on  either  side  of  the  focus;  a  similar 
arrangement  with  the  two  boundaries  assymetrical;  and  a 
circular  area  of  higher  velocity  to  one  side  of  the  fault. 
A  complete  set  of  records  for  both  P  and  S  waves  has  been 
obtained  for  the  six  different  models  for  each  of  the  four 
source  systems  under  study.  This  is  by  no  means  all  the 
cases  which  could  be  studied  profitably,  but  it  does  serve 
to  illustrate  the  type  of  result  to  be  expected. 
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Evaluation  of  the  P-wave  data  is  well  under  way. 

To  date,  the  shear-wave  records  have  been  studied  only 
for  the  homogeneous  model.  Table  I  summarizes  the  present 
status  of  the  analysis  of  the  records.  An  evaluation  of 
this  data  will  complete  phase  I  of  the  proposed  research 
and  will  be  prepared  as  a  report  for  presentation  before 
a  professional  society  and  for  publication.  Further 
experiments  may  be  directed  at  specific  special  cases 
differing  from  the  general  cases  covered  here. 

Preliminary  Results; 

In  general,  the  results  in  the  uniform  models  are 
consistent  with  theory.  The  amplitudes  and  directions  of 
first  motion  for  the  congressional  waves  for  all  uniform 
models  and  for  the  shear  waves  for  two  uniform  models  were 
presented  in  earlier  reports.  The  basic  patterns  for  the 
shear  waves  for  the  double-couple  and  radial-force  models  in 
a  uniform  disc  are  shown  in  Fig.  1. 

The  strength  of  the  shear  wave  was  found  to  be  greater 
than  that  of  the  compressional  wave  for  all  the  uniform 
models  except  In  the  case  of  the  radial  force.  For  the 
single  force  and  double  couple,  the  maximum  amplitude  of 
S  Is  about  two  times  as  strong  as  the  maximum  P.  For  the 
single  couple,  maximum  S  Is  about  four  times  as  strong  as 
p.  In  the  case  of  the  radial  force,  the  amplitude  of  S  Is 
very  small  compared  to  the  average  value  for  P.  Also,  from 
setting  up  numerous  radial-force  models.  It  appears  that 
the  smaller  the  scatter  In  P-wave  amplitudes  with  azimuth 
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the  smaller  the  S-wave  amplitudes  for  the  uniform  model. 
There  appears  to  be  a  definite  correlation  between  the 
overall  assymmetry  In  the  P-wave  amplitudes  and  the 
generation  of  the  S- waves. 

The  remainder  of  the  S-wave  data  will  be  discussed 
at  a  later  date.  Difficulties  have  been  encountered  In 
the  Interpretation  of  the  S-wave  first  motions  for  the 
non-uniform  models  due  to  the  presence  of  energy  arriving 
Just  before  and  during  the  arrival  of  the  main  shear  pulse. 
This  problem  is  also  encountered  in  the  interpretation  of 
actual  earthquake  seismograms. 

The  peak  amplitudes  for  the  compressional  waves 
in  the  four  uniform  models  were  measured.  A  consistency 
was  found  to  exist  between  the  first,  second  and  third 
peaks.  This  is  illustrated  for  the  case  of  the  single 
force  in  Pig.  2.  In  the  upper  graph  the  relative  values 
of  the  first  four  peaks  are  plotted,  a  positive  value  indicat¬ 
ing  up  on  the  record  (away  from  source  on  model).  The 
amplitudes  of  the  fourth  peaks  are  not  consistent  with  the 
second  peaks  due  in  part  to  the  simultaneous  arrival  of 
a  new  phase.  The  lower  graph  shows  the  normalized  peak- 
to-peak  value  for  the  first  two  peaks  (P^)  and  the 
second  and  third  peaks (P2_3).  The  S-wave  data  will  be 
studied  in  a  similar  fashion  as  will  all  the  data  for 
the  non-uniform  models. 
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The  presence  of  the  various  velocity  discontinuities 
have  not,  except  in  one  case,  reversed  the  direction  of 
first  motion  of  the  P  wave.  This  is  illustrated  for  three 
models  for  the  case  of  the  single-couple  and  radial-force 
systems  in  Pigs.  3  and  4,  respectively.  The  solid  curves 
in  Pig.  3  are  the  theoretical  curves  for  the  P-wave  distribu¬ 
tion  for  the  uniform  single-couple  model. 

Reversals  of  the  first  motion  of  P  compared  to  the 
theoretical  direction  for  the  uniform  model  are  present  in 
one  model  over  a  limited  azimuthal  range  at  the  boundary 
of  the  velocity  discontinuity  (Pig.  5).  This  model  had 
a  double-couple  source  with  a  velocity  discontinuity  running 
through  the  source  parallel  to  a  shear  plane.  The  reversals 
appear  to  be  due  to  the  initial  arrival  of  a  refracted  wave 
of  opposite  sign  from  that  of  the  direct  wave  at  points 
within  30  degrees  of  the  boundary  in  the  low-velocity  zone. 

In  the  cases  of  the  single-  aid  double-couple  models 
there  appears  to  be  a  shifting  of  the  nodal  points  in  the 
low-velocity  channel  of  less  than  15  degrees  toward  one 
boundary  of  the  high-velocity  medium. 

The  effect  of  a  boundary  of  any  kind  is  to  reduce 
the  amplitudes  of  the  P  waves  received  in  the  high-velocity 
medium  and  in  some  cases  in  the  shadow  of  the  high-velocity 
zone.  There  is  a  sudden  decrease  in  amplitude  upon  crossing 
the  velocity  discontinuity  into  the  high-velocity  medium  as 
can  be  seen  in  Pigs.  3  and  4.  Records  obtained  on  the 
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boundary  may  exhibit  normal  or  reduced  amplitudes  depending 
on  whether  the  pick-up  transducer  overlaps  more  of  the  low-- or 
of  the  high-velocity  medium.  A  quantitative  analysis  of 
the  amount  of  reduction  in  amplitudes  to  be  expected  due 
to  the  presence  of  the  boundary  and  the  difference  in 
elastic  properties  of  the  two  media,  will  be  made  later. 

II.  actual-displacement  experiments 

Frozen-fault: 

Figure  6  shows  a  graph  of  some  relative  amplitudes 
of  the  P  wave  for  the  actual -displacement  model.  The 
records  corresponding  to  these  data  were  obtained  in  the 
manner  described  on  page  9  of  Progress  Report  No.  3.  The 
azimuths  falling  within  the  shaded  areas  correspond  to 
azimuths  which  are  shielded  from  direct  radiation  from 
the  source  by  slits  (refer  to  to  Fig.  2  of  Progress  Report 
No.  3  for  geometry  of  the  source  region).  The  cosine 
curve  superposed  on  the  data  is  that  of  the  theoretical 
relative  amplitudes  to  be  expected  from  a  single-force 
source.  The  data  for  the  azimuths  0-180°  might  be 
expected  to  follow  this  theoretical  curve  more  closely 
than  any  other  because  the  slits  block  radiation  from 
the  tongue  side  of  the  model. 

Examination  of  the  graph  will  reveal  that  the  data 
do  not  follow  this  theoretical  curve  closely.  Although 
dilatation  might  be  expected  in  the  90-180°  region, 
compressions  are  observed  at  all  azimuths.  The  reason 
for  this  has  not  been  determined. 
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Mechanically  displaced  fault; 

The  results  over  the  past  months  obtained  using  the 
frozen-fault  technique  have  not  been  entirely  satisfactory 
and  the  method  is  very  time  consuming,  in  an  effort  to 
accellerate  the  accummulation  of  data,  other  techniques 
have  been  attempted  which  do  not  require  freezing.  One 
technique  which  shows  considerable  promise  employs  the 
same  type  of  model  and  electrical  circuitry  as  has  been 
heretofore  used,  but  stress  is  applied  simply  by  inserting 
a  blade  into  one  of  the  slits  of  the  model  and  twisting  it 
manually.  The  fault  surfaces,  being  in  tight  contact 
resist  slipping  because  of  static  friction.  Stress  builds 
up  on  both  sides  of  the  fault.  When  the  static  friction 
is  finally  overcome,  the  fault  surfaces  move  suddenly,  thus 
releasing  part  of  the  stress  in  the  form  of  elastic  waves. 

Records  of  the  radial  component  at  5  azimuths  made 
by  the  new  technique  are  shown  in  Fig.  J.  There  is  a 
clear  similarity  between  these  records  and  those  made 
by  the  frozen-fault  technique  (refer  to  Fig.  11,  Progress 
Report  No.  3).  The  P-wave  first  motion  appears  as  a 
distinct  pulse  in  both  cases.  There  is  also  a  strong 
long-period  component  associated  with  the  records  made 
using  both  techniques.  Reproducibility  of  records  is 
shown  for  two  azimuths  in  Fig.  7. 

Figure  8A  shows  examples  of  "double  breaks",  a 
phenomenon  that  is  also  prevalent  in  the  frozen-fault 
records.  In  this  case  the  stress  is  released  in  a 
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series  of  steps  resulting  in  overlapping  wave  trains 
(see  discussion  on  page  3  of  Quarterly  Report  No.  2, 
fourth  paragraph). 

The  records  of  the  two  methods  also  differ  in  some 
respects.  Although  the  wave  trains  appear  the  same,  the 
plots  of  direction  of  first-motion  against  azimuth  do  not 
appear  to  follow  the  same  pattern.  P-wave  motions  for  the 
frozen-fault  were  compressed  at  all  azimuths,  whereas  the 
new  technique  yields  compressions  in  the  range  0-90°  and 
dilation  from  90~l80°.  The  latter  distribution  is  con¬ 
sistent  with  the  hypothesis  that  the  source  mechanism  is 
equivalent  to  a  single  force. 

Radial  component  records  using  the  new  technique  have 
been  made  at  15°  intervals  around  the  model.  First-motion 
relative  amplitudes  have  been  plotted  in  Figure  9.  The 
theoretical  curve  for  the  single-force  source  has  been 
superposed  on  this  plot.  For  azimuths  0-180°  the  directions 
of  first  motion  agree  with  the  hypothetical  curve,  but 
amplitudes  are  widely  scattered  about  the  curve.  For 
azimuths  I8O-36O0  the  directions  of  first  motion  do  not 
always  agree  with  the  hypothetical  curve.  This  is  not 
surprising  since  most  of  the  recording  azimuths  in  this 
half  of  the  model  are  shielded  by  the  slits  from  the  direct 
P  wave.  Also,  the  presence  of  the  slits  will  cause  stress 
concentrations  which  could  be  expected  to  alter  the  first- 
motion  pattern  appreciably.  Another  difficulty  which  is 
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encountered  az  many  azimuths  from  18O-3600  is  the  first 
motion  is  not  a  clean  pulse,  but  rather  appears  to  be  made 
up  of  a  multiplicity  of  pulses.  The  interference  of  these 
pulses  makes  interpretation  of  amplitudes  at  these  azimuths 
ambiguous.  A  sample  record  demonstrating  this  interference 
effect  is  shown  in  Pig.  8B. 


TABLE  I 

Present  Status  of  Evaluation  of  Simulated  Model 

Data 


Single 

Force 

Single 

Couple 

Double 

Couple 

1.  Uniform  Model 

(1st  P-wave  peak) 

P,  P 

P 

p 

(1st  shear  peak) 

sf  P 

P 

N 

(2nd  and  3rd  P) 

Pp,Pc,  N 

R 

R 

2.  1/2  Aluminum  model 

(1st  P)  p  R 
(1st  shear)  S^  P 
(2nd  and  3rd  P)  ?2’?3  R 


P 

P 

F 


R 

F 

P 


3«  3/8  Aluminum  model 

(1st  P)  p  R 

(1st  shear)  S^  F 


N  R 

F  F 


4-*  3/2  low-velocity  zone 


(1st  P) 

Pl  R 

R 

R 

(1st  shear) 

sf  F 

F 

F 

(2nd  and  3rd  P) 

P2'P3  R 

R 

F 

5. 

3/3  low-velocity 

zone 

(1st  P) 

P.  R 

N 

R 

(1st  shear) 

Sf  F 

F 

F 

(2nd  and  3rd  P) 

P£,P3  R 

F 

F 

6. 

Spot 

(1st  P) 

P-,  R 

N 

R 

(lst  shear) 

F 

F 

F 

key:  P  =  in  a  previous  report 

N  =  in  this  report 
R  =  studied  but  not  presented 
F  =  records  made  but  not  interpreted 


Radial 

Force 


P 

N 

R 


P 

F 

F 


N 

F 


R 

F 

F 


N 

F 

F 


N 

F 


Fif*  I®.  Observed  amplitudes  of  radial  (P,  outer  curve)  and  transverse  (S,  inner 
curve,  ♦  indicates  clockwise  motian)  notions  far  a  radial-farce  source. 


jjlgrwga 


•mmmm ■■■■■■ 


■■■ ■ « w m i m m ---- - - ZZZZ --- 


Tig*  3*  Observed  P*erare  amplitudes  far  the  three  nonnudfapn  wodels  shewn, 

single -couple  source.  Carrs spooling  theoretical  curre  far  the  unlfors 
■edel  Is  shown  far  comparison.  Shaded  arias  represent  the  high-velocity 
sections* 


Fig.  U.  Obeerred  P-nre  aMplitudea  for  the  three  non-unifarB  modele  shown, 

Tadlul-farce  eouroe.  Shaded  areas  represent  the  high-velocity  seotione* 
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•  5*  Observed  P-vave  amplitudes  in  the  presence  of  a  velocity 
discontinuity)  double-couple  source.  Theoretical  curve  for 
uniform  model  is  shewn  for  comparison.  Shaded  area  represents 
the  high-velocity  section. 


Fig.  8a:  Double  breaks 


Enlarged  view 

Fig.  8b:  First  break  consisting  of  two 
overlapping  pulses 

Fig.  8:  Records  of  compound  pulses 
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